Selection of defect structures in twist-grain-boundary-A phase of chiral liquid crystals 
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We study the structure of the twist grain boundary in chiral smectic liquid crystals using the 
Landau-de Gennes model. By considering spatial variation of the smectic order, we distinguish the 
Melted-Grain-Boundary-A (MGBA) structure with cholesteric-like domains and the Twist-Grain- 
Boundary-A (TGBA) structure consisting of screw dislocations. The MGBA structure becomes 
stable near the transition to the cholesteric phase. On approaching the transition, the cholesteric- 
like domain grows outside the grain boundary. Also, the dislocation spacings of the TGBA structure 
agree better with experimental results than previous theories. 

PACS numbers: Valid PACS appear here 



Nontrivial physics in condensed matter is often caused 
by resolution of frustration. For example, mismatch be- 
tween the triangular crystalline order and antiferromag- 
netic spin order disturbs an equilibrium ordered phase 
even at the absolute zero [l[ ■ In the glass system compris- 
ing two types of hard core particles with different radii, 
the positional order of the particles and the bond order 
defined by the number of the nearest neighbor particles 
are frustrated, and generate disclinations [2|. 

Liquid crystals exhibit a wide variety of frustration- 
induced defect structures 0, B| • One of the reasons is the 
various kinds of order they show: nematic order, smec- 
tic layering order, and cholesteric (helical) order due to 
molecular chirality. Among them, frustration between 
the smectic and cholesteric order is paid great atten- 
tion because of scientific interest and potential applica- 
tion. Twist-grain-boundary (TGB) phase is the simplest 
one-dimensional defect phase, comprising series of finite- 
length smectic slabs (grains) with the layer normals ro- 
tating at a constant twist angle a, intervened by grain 
boundaries (Figflja)). It was predicted in the pioneering 
work [j| in terms of the analogy between liquid crystals 
and superconductors, and confirmed by experiment Q. 
As temperature increases, the ordered lattice of the screw 
dislocations "melts" into the dislocation fluid, called the 
chiral line (NL*) phase At higher temperature, de- 
fects construct a variety of three-dimensional networks, 
called smectic blue phases Q. Although the phase tran- 
sitions between these defect phases are observed, their 
spatial structures are not yet well understood. 

The grain boundary in the TGBA phase, having the 
molecular orientation parallel to the layer normal in the 
smectic slab, has been thought to consist of parallel screw 
dislocations, in analogy to the magnetic fluxes in the vor- 
tex lattice phase of superconductors. This structure (de- 
noted by the TGBA structure) was confirmed with the 
transmitted electron microscopy (TEM) observation [|. 
The length of the smectic slabs 4 and the interval of the 
screw dislocation Id were experimentally estimated [To| . 
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FIG. 1: Schematic representations of the TGBA phases, (a) 
The TGBA structure consists of the parallel screw disloca- 
tions (bold lines) . (b) The MGBA structure contains no screw 
dislocation. 



and theoretically calculated [ll|, ll2(]. However, the theo- 
retical values are about 20 times higher than the exper- 
imental values. This discrepancy in the defect spacings 
might be because of the theoretical treatment of the de- 
fect energy. 

Recently, another possible structure of the grain 
boundary has been proposed [HI]. In the TEM study, 
the screw dislocation array is not observed, which im- 
plies that the smectic order (and hence the disloca- 
tion array) is melted in the grain boundary (FigfT^b)). 
This structure is called the melted-grain-boundary-A 
(MGBA) structure, and is distinguished from the TGBA 
which has the periodic structure in the grain boundary. 

In this Letter, we consider the spatial variation of the 
smectic order, to analyze the structure of and compare 
the free energy of the TGBA and MGBA states. As 
by-products, we will calculate the characteristic lengths 
lb and Id in the TGBA state, and fix the discrepancy 
between theory and experiment. 

We consider one pair of the adjacent grain boundaries, 
and the smectic slab between them (FigfTJ). Setting the 
x-axis along the twist axis .we assume that the director 
n lies only in the yz-plane [5| , and express it by the local 
twist angle uu measured form the middle of the two grain 
boundaries as n = (0, — sin (uj + a/2) , cos (ui + a/2)). 



2 





213--., 


ot/2.-;- 





1® 



as 



cos(a/2) 



FIG. 2: The layer geometry of the TGBA structure seen from 
the twist axis (x-axis). The screw dislocations (bold lines) are 
supported by the smectic layers in front of the grain bound- 
ary (dot-dash lines) and behind the grain boundary (dashed 
lines). 



The smectic order is expressed by a complex order param- 
eter "J, and the equilibrium wave number of the smectic 
density wave is denoted by go- Layer displacement u is 
represented as ^ = ip exp [iqg (z — u)], where the z-axis is 
set along the layer normal at x = 0, and ip is the strength 
of the smectic order. 

The equilibrium structure is obtained by minimization 
of the Landau-de Gennes free energy consisting of the 
three components: 
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The smectic energy F am adjusts |^| to the equilibrium 
value ipo- Ffi is the Frank elastic energy of n, where 
the three elastic moduli of the splay, twist, and bend 
components are set to be equal, fco is the inversed helical 
pitch due to the molecular chirality. The coupling energy 
F cp is the cross term of ^ and n. The correlation lengths 
of ^ (n) is £ (A), proportional to 1 / v^NA — T. Tna is 
the nematic (N)-smectic A (SmA) transition temperature 
for the achiral case fc = 0. 

For the TGBA, the smectic order is destructed near 
the dislocations by the helical order. In the grain bound- 
ary with the thickness ~ £, the director changes slowly 
so that n ~ e z (FigEJ, because £ <C A in typical experi- 
ments 13]. The layer displacement is approximated with 
the linear superposition of the screw dislocations [l4[ 
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This layer distortion characterize the screw dislocation 
array, and F cp , containing u, is dominant near the grain 
boundary. Close to the dislocation cores, the local profile 
of "0 is a solution of 8F cp /Sip = 0, which is approximated 



v *=u 

The solution is 
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ip = C ^ cosh (2nx/ Id) — cos (27ry/Z d ), 



if,. (5) 
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where C is an integration constant [Tja]. Near the SmA- 
TGBA transition, ip increases up to ipo outside the grain 
boundary x > £. The energy density cost of making dis- 
location is e = 27ri/>Q^a/2, and the energy density gain 
by twisting the layers is ipgqQA 2 koa. Equating the two 
contributions, we determine the lower critical field for the 
SmA-TGBA phase transition ko = fc cl = l/2qoA 2 , which 
is much lower than the previous Renn-Lubensky result 
kf^ = ln(A/£) /2<7qA 2 [g]. In our theory, the energy per 
dislocation is lowered, because the long range layer dis- 
placement around the single dislocation (|Vu| cx 1/r) is 
summed out (Q}, and also because of the relaxation of 
ip whose correlation length £ is shorter than that of the 
director A, which has been neglected in previous theories. 

For the MGBA structures, ip is uniformly small in the 
grain boundary. The solution of the equilibrium condi- 
tion 5F cp /Si/j — is 



ip = D sin — , 



(7) 



where D is an integration constant [la ] . 

The energy of the grain boundary thus can be calcu- 
lated with © for the TGBA, and for the MGBA. It 
can be shown that the director rotates at a constant rate 
because of the weak smectic order. 

For both the TGBA and MGBA, the smectic or- 
der grows and F sm becomes one of the main contribu- 
tions at distance larger than £ from the grain bound- 
ary. In the smectic slab x > £, the twist is expelled, 
and the layer normal e z + Vu is close to a constant 
(0, — sin a/2, cos a/2). Using this relation for u, the sum 
of the smectic and coupling energies is reduced to the 
effective smectic energy 
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where ip is the effective equilibrium smectic order, and 
u> ranges from —a/2 to a/2. Because of the frustra- 
tion between the smectic and helical order, ip is reduced 
at large u>. According to ©, if w increases beyond 
a c /2 = 2 arcsin (1/2<7oC)j the equilibrium smectic order 
vanishes, and a cholesteric-like domain appears near the 
grain boundary. In the cholesteric domain, the director 
rotates at a constant rate, which we denote u)' . 

If the twist angle a is smaller than a c , ip approxi- 
mately equals ip — ip(uj) in the smectic slab x > £, 



3 



because the characteristic length scale of u> and /p is A 
(3> £). The twist of the director penetrates the smec- 
tic slab with the characteristic length A in the form 
uj oc sinh [(x — i&/2) / A]. Combining the solutions for the 
grain boundary and the smectic slab, we obtain the total 
free energy densities for the TGBA and MGBA states, 



/tgba/V'o 
and 
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where the dimensionless chiralities Pt — (k — k c \)/k c 
and Pm = k/k c are introduced, and X = (— £ + 4/2) /A 
is the ratio between the length of the smectic domain and 
the director correlation length. 

Minimization of the free energy density (fTP|) deter- 
mines lb and Id = d/2sm(a/2) for the TGBA state 
as a function of the thermodynamic variables T and 
kg, To the lowest order of K , the equilibrium struc- 
tural lengths are obtained as lb/ X = 2(X + k _1 ), 
and ld/X — ^/2tt/ Pt (1 + Ktanh X), where kX = 



1 + V/W3/3+ {/I-VI-/ 3 / 3 



and /3 



(Pt/k) 2 - In typical experiments d, OH, -Pt ~ 0.1, and 
k ~ 100. The ratio remains on the order of 1, 

agreeing with both the experimental and previous theo- 
retical results [13, El, E3 ■ The order of lb/\ and Id/ A is 
0.1 — 0.01 except near the SmA-TGBA phase transition, 
which agrees again with the experimental results [Ic| . In 
our theory, introducing the spatial variation of ip, the 
defect energy is evaluated much lower (so that the equi- 
librium defect density is larger) than in the previous the- 
ories [n m. 

Change of the stability between the TGBA and MGBA 
states occurs when /tgba = /mgba- This condition 
is satisfied for Pt/k,Pm/k ~ 0.10, which corresponds 
to the twist angles a ~ 0.41a c for the TGBA and 
a ~ 0.56a c for the MGBA. Averaging the conditions, 
in experiment, a TGBA- like structure would change to a 
MGBA-like structure at a = c*tm — 0.5a c , or Id — 6£. 
This means that the MGBA can be understood as over- 
lapped screw dislocations each of which has core radius 
on the order of £. In a real situation, transition between 
the TGBA and MGBA might not be thermodynamic, 
and could be a crossover with a gradual change of the 
grain boundary structure. Note that here we compared 
the two limiting cases (the TGBA with the linearly su- 
perposed dislocations and the MGBA with the perfectly 
uniform grain boundary). 

If the twist angle is large such that a > a c , a local 
cholesteric-like domain appears where \lo\ > a c /2. Ex- 



istence of the cholesteric domain means that the smec- 
tic order melts at the whole grain boundary. Thus if 
a > a c , the structure is always MGBA. We obtain 
an approximated equation for tp in the cholesteric do- 
main as tp" — q^tp (^'qX ~ a /2) 2 — W£ 2 j with the help 
of sin(w/2) ~ J/2 for \u\/2 < tt/4, and g 2 £ 2 > 1 
for the typical case [Uj]. In the London limit 0, [Tl| . 
using A 3> £ ^> Qq 1 , the solution is approximated as 

if) oc |t7|[^ w ° 90 ^ ) 1 ] ^ 2 exp (— r) 2 /2) , where r\ is the re- 
duced coordinate (a; — a/2uj' Q ) /£ with the effective cor- 
relation length £ = 1/ y/qoujQ. Defining L by uj (x = L) = 
— a c /2, the cholesteric domain extends to x = L + £, be- 
cause Tp (L + £ + Sx) ~ tpo (qoS,) 1 ^ 2 [(£ + fix) /M 1 ^ 2 grows 
to the order of ipo as Sx increases to the order of A. On 
the other hand, the smectic domain satisfying ip = ip 

ends at x — L + £, because W (L + £) /tpo] ~ £ and 
the gradient term in ([5]) becomes a major contribution. 
Combining the three domains, the grain boundary, the 
cholesteric domain, and the smectic domain, we obtain 
the total free energy density. The smectic domain is van- 
ished at the transition to the cholesteric (N*) phase, and 
the upper critical field fc C 2 is shown to have the same value 
l/qo£, 2 as in the Renn-Lubensky result [f|. The length 
of the grain boundary and cholesteric domain is plotted 
in Figf5Ja). The cholesteric domain starts to grow more 
rapidly than A at the temperature where the twist angle 
a equals the critical angle a c . Thus if the twist angle for 
this domain growth is measured, one can indirectly esti- 
mates the characteristic angle of change from the TGBA 
to MGBA otm — 0.5a c . Note, the abrupt growth of 
the cholesteric domain in Fig[3]Ja) might be an artifact. 
Instead, the growth should be continuous [161 ]. 

The phase diagram is shown in Fig[3Jb). The equilib- 
rium state changes as SmA-TGBA-MGBA-N* with the 
increment of the temperature and chirality. This state 
sequence can be understood in terms of symmetry. Both 
of the SmA and TGBA have the discrete translational 
symmetry of the smectic order in the yz-plane. Only the 
pair of the TGBA and MGBA possesses the discrete he- 
lical symmetry along the x-axis. For the MGBA and N*, 
they are only the pair having the continuous symmetry in 
the yz-plane. The state pairs having the common sym- 
metry are neighboring in the phase diagram. In addition, 
stability of the TGBA state to the SmA phase is higher 
than that in the previous work, because of the energy 
cost of the grain boundary is lowered by the interaction 
between dislocations in each grain boundary, as well as 
by the relaxation of the smectic order parameter. Such 
changes of the defect energy and structure by introduc- 
ing the spatial variation of the smectic order may occur 
also in the TGB-embedded double twist cylinder [13] and 
other smectic blue phases [|[ . The TGB phase has a sim- 
ple one-dimensional defect structure, and our results may 
provide a fundamental clue to understand the structures 
of these complex defect phases. We also note that the 
MGBC* phase was shown to be always stable over the 
TGBC* phase due to the macroscopic helical order of the 
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director along the layer normal The MGBA struc- 
tures does not have such a helical order and its stability 
over the TGBA was far from trivial. 
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FIG. 3: (a) The temperature dependence of the size of the 
grain boundary and cholesteric domain (divided by A). The 
temperature range between the SmA-TGBA and MGBA-N* 
transition is shown. (b)The phase diagram. Solid lines mean 
the second order phase transition. The dot-dash line shows 
the previous result for the SmA-TGBA transition Dashed 
lines separate different equilibrium structures. 

Let us discuss our results in terms of the analogy be- 
tween the TGB phase and the vortex lattice phase of 
superconductors. The smectic order parameter ^> and 
the director n correspond to the wave function of the 
superconducting particles and the vector potential, re- 
spectively. The screw dislocation in the TGB phase, the 
twist angle a, and the molecular chirality ko correspond 
to the vortex, the magnetic flux, and the external mag- 
netic field. The difference is the dimensionality: In the 
TGB phase, three spatial dimensions are divided into 
one for the direction of the screw dislocation, one for the 
twist axis, and one for the periodic dislocation array. In 
the vortex lattice phase of isotropic superconductors, on 



the other hand, the three are divided into one for the 
direction of the vortex, and two for the lattice period- 
icity. However, in anisotropic superconductors, vortices 
are aligned in a one dimensional chain for a certain mag- 
netic field range as in the TGB phase. In liquid crystal, 
the origin of the corresponding anisotropy is the molec- 
ular chirality. Theoretical study of the vortex chain is 
based on the interaction of the phase of the wave func- 
tion, corresponding to the layer displacement in liquid 
crystals Thus it is interesting to introduce the spa- 
tial variation of wave function near the vortices, as in 
our theory for chiral liquid crystals. The vortices in a 
chain are "melted" to the Josephson vortex by stronger 
anisotropy. Thus one might regard the MGBA state as 
an analogue to the Josephson vortex state. 

In summary, we analyze the TGBA and MGBA states 
of chiral liquid crystals with the Landau-de Gcnncs 
model. By considering the spatial variation of the smectic 
order parameter, we found that a cholesteric-like domain 
appears at twist angle higher than the critical angle a c . 
We have shown that the proposed MGBA state [l3[ is cer- 
tainly stable over the TGBA state when the twist angle is 
higher than cvtm — 0.5a c . We may indirectly specify the 
MGBA state through observing the cholesteric domain 
that appears at a = a c . The energy cost of the dislo- 
cation is also greatly changed from the previous theory, 
resulting in an agreement of the dislocation spacings lb, 
Id with the experiments (hJ, and the enhanced stability 
of the TGBA state over the SmA phase. It would be 
interesting to extend our results to other defect phases, 
including the smectic blue phases of liquid crystals, and 
the vortex chains of anisotropic superconductors. 
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